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1.  INTRODUCTION 

This  document  discusses  a  first-order  vector  perturbation  theory  calculation  for 
angle-resolved  scattering  (ARS)  from  multilayered  optical  components  having  random 
roughness  at  each  interface  of  the  composite  stack.  The  theoretical  basis  of  this  work  was 
first  developed  and  presented  many  years  ago  (Reference  1).  Comparison  between  theory 
and  experiment  yielded  good  results  (Reference  2).  Since  then,  technological  devices 
involving  specialized  multilayer  optical  components  having  many  tens  or  hundreds  of  layers 
have  been  conceived  or  fabricated  that  depend  critically  on  very  high  levels  of  performance 
such  as  99.99%  reflection  or  transmission.  This  means  that  scattered  light,  which  is  on  the 
order  of  0.0 1  %  or  more  of  the  incident  energy  either  in  reflection  or  transmission,  may  limit 
or  destroy  intended  performance.  Thus,  knowledge  of  expected  scattering  characteristics 
for  multilayer  components  in  optical  systems  definitely  can  be  a  valuable  design  tool.  This 
theory  can  be  applied  also  to  uncoated  components. 

Also,  a  Fortran  source  code  for  numerical  evaluation  of  the  theory  presented  here  is 
given  along  with  user  instructions  and  sample  output  data.  Formulas  for  reflection  and 
transmission  scattering  are  given  and  are  related  to  the  source  code. 


2.  BASIC  THEORY 


The  multilayer  stack  is  assumed  to  be  illuminated  by  a  monochromatic  plane  wave  of 
wavelength  X  incident  at  polar  angle  Sq  relative  to  the  direction  normal  to  the  mean  surface. 

Regarding  polarization  of  the  incident  beam,  the  electric  vector  can  be  oriented  at  an  arbitrary 
angle  relative  to  the  plane  of  incidence.  Referring  to  Figure  1,  interface  j  of  the  multilayer 

stack  is  described  by  z  =  dj  +  AZj(p),  where  Az/p)  is  a  random  variable  that  describes  the 

surface  roughness  of  interface  7  and  fluctuates  about  the  mean  surface  level  dj.  The  roughness 

of  the  interface  j  has  a  mean  square  value  of  (|  Az^(p)  p)  =  6^  when  averaged  over  the  interface 

j,  and  p  =  (x,  y )  is  a  point  on  the  mean  surface  reference  plane  of  interface  j.  Since  the 
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FIGURE  1.  Schematic  Showing  Nomenclature  Associated  With  the  Multilayer  System  That  Consists  of  a 
Substrate,  L  Layers,  and  a  Superstrate.  The  superstrate  contains  the  incident  beam  specular  beam,  and  reflected 
scattered  field.  As  shown,  the  substrate  contains  only  the  transmitted  scattered  field,  if  any,  although  there 
may  also  be  a  transmitted  specular  beam.  The  angles  Sg  and  6  are  the  polar  angle  of  incidence  and  scattering, 
respectively.  The  incident  beam  is  assumed  to  be  in  the  (x,z)  plane.  The  scattered  field  has  direction  fG.tp) 
where  the  azimuth  angle  tp  is  not  shown  here.  The  regions  of  the  multilayer  system  have  complex  dielectric 
constants  Cy  where  j  =  1  is  the  substrate,  j  =  2 -*L  + 1  are  the  layCTS  of  the  stack  with  mean  physical 
thicknesses  Xj,  and  j  =  L +  2  is  the  superstrate.  The  dj  are  the  cumulative  thickness  of  the  multilayer  stack 
at  the  mean  surface  level  of  interface  j,  and  the  function  Az/z,  y)  is  the  random  roughness  function,  which 
deviates  about  z  =  dj. 
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results  given  in  this  work  are  based  on  a  first-order  perturbation  calculation,  it  is  assumed 
that  the  ratio  of  root-mean-square  (rms)  roughness  to  wavelength  bj/X  is  much  less  than 
unity. 

The  multilayer  stack  has  L  layers,  which  are  labeled  j  =  2  -¥L+  1.  The  total  multi¬ 
layer  system  is  assumed  to  consist  of  these  L  layers,  each  having  physical  thickness  Zj  plus 

a  semi-infinite  substrate  and  superstrate.  For  convenience  we  will  apply  the  term  "layer" 
to  the  substrate  and  superstrate,  where  the  substrate  is  layer  j  =  1  and  the  superstrate  is 
layer  j  =  L  +  2.  Thus,  the  total  multilayer  system  has  L  +  2  layers.  The  homogeneous 
complex  dielectric  constant  of  layer  j  is  Ej,  where  ;  =  i  L  +  2.  In  this  work,  the  real  part 
of  Ej  can  be  positive  or  negative,  but  the  imaginary  part  of  Ej  is  always  positive. 

In  each  layer  of  the  multilayer  system,  we  seek  electric  and  magnetic  field  vector 
solutions  of  the  wave  equation  having  tangential  components  that  are  continuous  to  first- 
order  in  bj/X  across  the  rough  interface  j.  As  will  be  seen,  this  is  equivalent  to  a  distribution 
of  surface  currents  at  the  planar  mean  interface.  For  interface  J,  the  unit  normal  to  first-order 
in  bj/X  may  be  written 


8Az/p)'|_/aAr/p)" 

dx  J  t  ay  , 


(2.1) 


In  order  to  satisfy  the  conditions  of  continuity,  the  electric  field  in  medium  j  may  be  written 
as  a  sum  of  zero-  and  first-order  fields,  E ^(p,  z) = ^”\p,  z ) + ^*\p,  z ).  This  may  be  expanded 

to  first-order  about  the  z  =  dj  plane  to  approximate  the  electric  field  at  z  =  dj  +  Az(p),  and 
this  yields 

%ip,dj  +  Mp))  =  '^\M)  +  ^jip)^'^\p,z)\.,,^^^  (2.2) 


■“^0) 

The  Ej  is  the  zero-order  field  solution  for  layer  j  that  results  when  all  interfaces  are  perfectly 

smooth,  i.  e.,  Az/p)  =  0  for  all  j.  The  Ej  is  the  first-order  correction  (scattered  field)  to  the 

total  electric  field  for  layer  j,  which  arises  as  a  result  of  interface  roughness.  Considering 
first  the  boundary  condition  that  the  tangential  components  of  the  electric  field  given  in 
Equation  2.2  be  continuous  across  interfaces  j  and ;  +  7,  we  have 

nxMj^ij(p,dj  +  AZj(p))  =  0  (2.3) 
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where 


LEj ,  dj  +  Az/p))  =  Ej ,  ,(p,  dj  +  Az/p))  -  Ejip,  dj  +  Az/p))  (2.4) 

is  the  difference  of  the  electric  fields  evaluated  across  the  z  =  dj  +  Az^(p)  boundary.  Using 
Equations  2.1  and  2.2  in  Equation  2.3  yields 


for  discontinuities  of  thex-  and  y-components  of^’^and^+j  across  interface  j. 


We  assume  that  the  zero-order  field  in  layer  j  has  the  form 
^®\p,  z)  =  z)e‘*“  ’  ^  =  TPikfi,  z)e‘^ 


(2.6fl) 


which  indicates  that  the  incident  beam  is  in  the  (x,z)  plane  where  the  vector 
ko  =  \ko  =  x(co/c)VeL+2sin  6o  is  the  projection  of  the  incident  wave  vector  onto  the  mean 
surface  plane.  Note  that  the  incident  medium  (j  =  L  + 2)  must  be  lossless. 


The  first-order  field  in  layer  j  is  written  as  a  Fourier  integral 


(2.6ft) 


where  k  =  {K,ky)  is  the  projection  of  the  scattered  field  wave  vector  on  the  (x,y)  plane.  An 

exp(-zcor)  time  dependence  is  assumed  but  not  shown  explicitly.  Using  Equations  2.6  in 
Equations  2.5  yields 


{lnf^efl,J^^{k ,  dy)  =  dj)^j{k,  -  k) 


(2.7a) 
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=  (2.76) 

where 

Az/r„-r)=JdV4z,(p)e'^''’'"  (2-8) 


These  equations  indicate  that  the  jc-  andy-components  of  the  first-order  electric  field  solution 
are  discontinuous  across  the  mean  planar  surface  at  z  =  dj.  This  implies  that  there  is  a  sheet 
of  Dirac  6— function  surface  currents  at  the  plane  z  =  dj  that  is  equivalent  to  the  roughness 
(Reference  3).  Equation  2.8  is  the  Fourier  transform  of  the  profile  for  interface  j. 

Similar  continuity  equations  for  the  x-  and  y-components  of  the  magnetic  field  yield 


and 


(2.9a) 


(2.9b) 


where  the  displacement  vector  =tjEj  and 


(2.10) 


The  first-order  scattered  field  is  obtained  primarily  by  finding  solutions  to  Maxwell’s 
equations  that  have  discontinuities  as  given  in  Equations  2.7  and  2.9.  This  requires  solution 
of  the  zero-order  field. 


3.  ZERO-ORDER  FIELD 


The  wave  equation  for  the  zero-order  electric  field  in  layer  j  of  the  multilayer  system 
is  written  as 

[^^  +  (G)/c)%]Ef\p,z)  =  0  (3.1) 

Assume  layer  L  +  i  is  illuminated  by  a  linearly  polarized  plane  wave  with  wave  vector  in 
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the  x-z  plane  incident  at  polar  angle  Gq  relative  to  the  z-axis.  The  incident  beam  may  be 

written  as  the  sum  of  a  p-polarized  (electric  vector  parallel  to  plane  of  incidence)  and 
s-polarized  (electric  vector  perpendicular  to  plane  of  incidence)  component.  Thus,  the 
zero-order  electric  field  amplitude  solution  to  Maxwell  ’  s  equations  for  layer y  has  the  general 
form 


^^p,z)=:e;°^(p,z)-h^;(5,2) 


(3.2) 


where  the  p-  and  s-polarized  component  solutions  are 


^°J,(p, z)  =  {xefliz)  -I- zef^liz)}  e 


(0)/ 


=  cosma]  ’ — - ‘  +b)' — - p=-c  '  >e 

[  (a)/c)^|EJ  (to/c)V^  J 


(3.3fl) 


and 


^!l(P>z)  =  -yej®;(z)e“**  =  -y sin\j;| 


itj* 


(3.3h) 


respectively.  The  angle  of  the  electric  vector  relative  to  the  plane  of  incidence  is  y.  Thus, 
y  =  0  degrees  is  p-polarized  incidence  and  y  =  90  degrees  is  s-polarized  incidence.  In  layer 
j,  we  have  {(O/cfej  =  k^  +  qf^  or 


qf^  =  {(a>lc)hj-kif 


(3.4) 


with  co/c  =  Infk.  To  compute  the  unknown  coefficients  af*,  bj°\  gf  \  and  jf\  we  equate 

tangential  components  of  the  electric  and  magnetic  field  amplitudes  across  the  interface 
boundaries,  use  the  facts  that  there  are  no  incoming  beams  from  the  substrate  region  and 
that  the  incident  intensity  is  of  unit  magnitude.  This  yields  solutions  of  the  form 


^(0)  p(l.>-l) 
(0)_  HL  +  V  0,12 


(3.5a) 


^(0)  p(l,y-l) 


0.5b) 
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^(.0)  o(l,7-l) 

(0)  Vi +2^^  0,12 

'  «/”«*« 

(3.5c) 

(0)  o(l,7-l) 

^0)  Vi+2^0,11 

(3.5d) 

h  ^(0)5(l.i.l) 

where  ;  =  1  — >L+2. 

For  j  =  2  — >  L  +  2,  the  matrix  elements  Po.Vi 

■  and  Pa,\2  are  derived  from  the  matrix 

product 


pO-J-i)- 


p  (1.7-1) 
~  0,11 

-1) 


>(1.7 -1)> 


p(1.7 
^  0,21 


0,12 

p(X,i- 

^0,22 


1) 


-Po 


(i-i) 


(3.6a) 


and  when  j  =  l,wc  define 


p(1.0)_ 
^0  ~ 


p(1.0) 

'0,11 

p(l.O) 

^0,21 


j(l.O)^ 

0,12 

jd.O) 

0,22  y 


ri  0) 

lo  ij 


(3.6b) 


For  m  =  l^L  +  l,the2x2  matrices  in  the  product  on  the  far  right  of  Equation  3.6fl  have 
the  form 


Po,!! 

\P0,2l 


Po,^J 


(3.1a) 


where  the  elements  are  given  by 


(m) 

Po,ll  = 


= 

P0,12 


2^iNemem+i 


exp[j(^ri  1-041 


expHCO  1  +  041 


(0)x 


(3.1b) 


(3.1c) 
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fp  4,(0)  _p  47(0) 
(,n)  ^ynHm  +  1  'hn  +  lHm 

^0.21  I  2<7®^e„e„+i 


expLi’C^i^l  1+041 


(3.7cf) 


fp  47(0)  4.0  47(0)^ 

pS=  ■—  expH(,“.-?rVj 

I  2<7„'V£m£m  +  l  . 


(3.7e) 


The  corresponding  s-polarized  field  matrix  elements  and  for  j  =  2—^L+2 
are  given  by  the  matrix  product 


rc(ij-i)  c(U-m 
0.11  _  (1)  (2)  0-1) 

*^0  -  ^(1,7-1)  ^(1.7-1)  “^0  -io  -"io 

V'jQ.Sl  *^0,22  J 


(3.8a) 


and  for  7  =  i,  we  define 


f  -(1.0)  -d.O) 

r(1.0)_  *^0.11  *^0,12 

*^0  -(1.0)  -(1.0) 

W0.2I  *^0.22 


1  0 
0  1 


(3.8b) 


For  m  =  l— >L  +  l,the2x2  matrices  in  the  product  on  the  far  right  of  Equation  3.8o  have 
the  form 


-V-/ 

■>0 


-(<«) 

-(«) 

•*0.11 

•*0.12 

-(«) 

-(«) 

•*0.21 

^0.22 

(3.9a) 


where  the  elements  are  given  by 


•^0.11  ~ 


^  exp[z(0 1-041 


(3.9b) 


/'4,(0)_-(0)  ^ 

^?12=  ■  expH(Oi+^«V«l 

I 


0.9c) 


f47(0)_-(0)  A 

4%=  ■  exp[i(0i  +  O41 

V  ^Qm  j 


0.9d) 
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-W  _ 
■»0,22~ 


V 


(0) 

m 


exp[-i(^rii-?rv*] 


(3.9e) 


In  Equations  3.5  for 7  =  L  +  2,  the  =Jl12  =  1.  which  validates  the  assumption  that 

with  Equations  3.3  the  amplitude  coefficient  for  the  incident  beam  in  the  superstate  is  of 
unit  magnitude.  Further,  with  Equations  3.5,  we  see  from  Equations  3.6b  and  3.86  that 
<2^ =  0  which  validates  the  assumption  that  there  are  no  incident  beams  from  the 
substrate  medium. 

At  this  point,  the  zero-order  field  is  solved  and  the  intensity  reflection  R  and  trans¬ 
mission  T  for  a  p-  or  s-polarized  incident  beam  are  given  by 


(3.10a) 

7',=ii>n'9!tor')«(9™2) 

(3.106) 

(3.10c) 

(3.10d) 

where  91(Q)  denotes  the  real  part  of  Q. 

The  average  flow  of  power  from  the  superstrate  region  incident  upon  the  j  =  L+  1 
interface  may  be  calculated  from 

(3.11) 

where  only  the  incoming  parts  of  and  are  retained,  i.  e.,  those  terms  proportional 

to  64+2  and  ^+2  ia  Equations  3.3.  The  integration  is  over  the  area  A  of  illumination,  which 
has  dimensions  assumed  to  be  essentially  infinite  compared  to  wavelength  X.  It  is  further 
assumed  in  this  work  that  the  dielectric  constant  e^, + 2  of  the  incident  medium  is  nonabsorbing, 

and  thus  qP+2  is  also  real. 
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4.  FIRST-ORDER  (SCATTERED)  FIELD 


Breaking  Equation  2.6b  into  p-  and  s-polarization  components,  the  first-order  solution 
for  the  scattered  electric  field  in  layer  j  is  rewritten: 


/  2  ^2 
ylTCj 


(4.1a) 


and  the  magnetic  field  is  obtained  from 


Vx^.‘\p,z) 


(4.16) 


where 


.‘•r.  .  (.,(fc£2iiy 


and 


?;;;(z)=(kxz){j(»e‘-,»>.^^.,^-<,; 


The  in-surface  wave  vector  of  the  scattered  field  is 


k=kk  =  ixcos^  +  y  sin  ^)k  =  (k^,  ky) 


(4.1c) 


(4.1d) 


(4.2a) 


where  0  is  the  azimuth  angle  of  scattering  relative  to  the  x  direction.  Since  we  require 
(co/cfEj  =  k^  +  qj^\  the  magnitude  of  the  z-component  of  the  scattered  field  wave  vector  is 


qf^=[iai/cfej-k^r  (42b) 

The  k  is  generally  different  for  scattering  into  different  media.  This  is  discussed  further  in 
the  following  section.  In  this  work,  we  only  calculate  scattering  into  the  reflection  (£^+2 
medium)  and  transmission  (Ej  medium)  hemispheres  rather  than  into  any  of  the  L  layers 
between  these  two  semi-infinite  media. 
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The  coefficients  5l+2  are  called  the  reflection  scattering  coefficients.  The 

and  are  the  transmission  scattering  coefficients.  These  coefficients  are  calculated 

by  application  of  Equations  2.7  and  2.9  with  Equations  4.1  and  can  be  written  after  some 
manipulation  as  a  sum  over  the  L  +  7  interfaces  as 


(1)  (m^c)  I:)  unuu-n 

Sl.i  2sjw..)^, 


1  I  +  l  -• 

A  M  ax 


[  (^^W^y(0)  (U-l) _ l_  y(0)  (l.y  - 1)  1, 

(1) 


f<i)  -  (^g)  y(0)r  <,(1.1, + i)e(i.7  - 1) _  o(i,i. + - 1)-. 

•fl  ~/,oa.L  +  nA  -(1)  92  *521  Si  / 


(4.3a) 


(4.3h) 


(4.4a) 


(4.46) 


The  ^j(ko  -  k)  is  the  Fourier  transform  of  the  interface  roughness  function  Az,(p)  as  given 

in  Equation  2.8.  The  Fi®],  and  Y^j  are  called  the  zero-order  source  terms  and  may  be 
written  as 


x;°^  =  iXoyc)(e^,i--ep{ej2^)sin<l)-e;>,.)cos(|)} 


YS  =  ik 


_p  ^ 

^j+i 


and 


yg  =  /(co/c)  (e,.,,  -  £,.)  {ef;iidj)  cos  <[)  -n  eflidj)  sin 


(4.5a) 

(4.56) 

(4.5c) 


To  evaluate  further  the  scattering  coefficients  in  Equations  4.4,  the  4,  T),  and  p.  functions 
are  required  and  these  are  defined  by 


Si  ~‘^l.ll  ^  ‘^1,12  ^ 


(4.6a) 
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S2  “■•^1,21  ^  *^1.22 

„a.y-i)  _  pd.y-D^'^i’H  .  pd.y-i)g^y^H 

Ml  1,11  ^  ^^1,12 


T,d.y-i) _  pd.y-i)."^i'H- _pd.y-iL‘«y^‘‘(r 
TI2  1.11  ^  ^1.12 


,d.y-i) _  pd-y-D^-’/H  + 


„ii.y-ii_  pu.y-‘v 
M3  ~^1.21 


1.22 


^d.y-i)^pd.y-i)^-wi’H_p(U-i)^w7H 


(4.6&) 

(4.7fl) 

(4.7fe) 

(4.7c) 

(4.7d) 

(4.8fl) 

(4.86) 


The  Ti  and  |i  functions  pertain  to  the  p-polarized  scattered  field  and  require  evaluation  of 
the  matrix  products 


d.y-i)_ 


( pdiy~i)  pd-y~i)') 
^1,11  ^1.12 

pd.y-i)  pd.y-1) 
V  1.21  “1.22  ) 


=pPV.®.-.rf-” 


(4.9a) 


where  ;  =  2  -4  L  +  2.  For  j  =  1,  we  define 

/^pd.O)  pd.oA  /  „N 
pd.O)^  M.12  ^1  0 

pd.O)  pd.O)  I  0  1  ) 


(4.96) 


where  the  2  x  2  matrices  in  the  product  on  the  far  right  of  Equation  4.9fl  have  elements 

|exp[i(ql'ii-ql*Vm3  (4.10a) 


= 

Pi.n 


= 

Pi, 12 


+1  J 


(p  -F 

^mHm  +  1  ^m  +  IHm 

2^ijVEmem  +  l 


expH(^mii  +  <7mVj 


(4.106) 
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= 

Pl,21 


exp[i(^i'ii+^i‘Vj 


(4.10c) 


expH(^i‘li 

J 


(4.10c?) 


For  the  ^  functions,  which  pertain  to  the  s-polarized  field,  the  corresponding  matrices  for 
y  =  2  -> L  +2  are  given  by 


V 


cUJ-iA 

•^i.n 

co.y-i) 

•Jl.22 


—  iJj  iJj 


(4.11a) 


For  y  =  1,  we  define 


fcd.O) 

‘*1,11 

cd.oA 

**l,12 

c(l,0) 

c(l.O) 

(,■*1.21 

•^1.22  y 

1  o' 
0  U 


where  the  elements  of  the  matrices  on  the  far  right  of  Equation  4. 1  la  are 


2«, 


(1) 


exp[i(^iJl, 


-(»«)  _ 
^1,12  ~ 


Hm  Hm^\ 


H 


(1) 


expH(^l‘^i  +  qiVm] 


-(«)  _ 
*1,21  ~ 


2q, 


(1) 


exp[i(^l‘li  +  ^l‘Vj 


^qm+qmh 


2q, 


(1) 

m 


lexpH(^l‘li-ql‘Vj 


(4.11&) 


(4.12a) 


(4.12h) 


(4.12c) 


(4.12d) 


These  and  Pf’-'  matrices  are  analogous  to  those  used  in  the  evaluation  of  the 

zero-order  field  coefficients  as  in  Section  3. 
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5.  ANGULAR  DISTRIBUTION  OF  SCATTERED  LIGHT 


We  seek  a  ratio  dPIdO.  of  power  dP  scattered  into  direction  (6,(j))  divided  by  solid 
angle  dQ  =  sin  QdQd^.  This  ratio,  normalized  by  the  total  power  incident  on  the  surface  Pq. 

is  the  ARS,  {\IP^{dPldO).  This  expression  is  related  to  the  bidirectional  reflectance 
distribution  function  (BRDF)  by 

ARS=BRDFcose  =  -^^  (5.1) 

PodQ  ’ 


Calculation  of  the  ARS  is  accomplished  by  considering  an  integral  of  the  form 


(5.2) 


where  j  =  L  + 2  for  reflectance  scattering  and  j  =  1  for  transmission  scattering.  For 
j  =  L  + 2  and;  =  i,  is  the  time-averaged  energy  propagating  away  from  the  multilayer 

stack  in  the  positive  and  negative  z  direction,  respectively.  Using  Equations  4.1  for;  =  L  +  2 
in  Equation  5.2  yields,  after  some  manipulation. 


(5.3a) 


and  likewise  for  j  =  1  we  obtain 


(5  3b) 


where  the  minus  sign  arises  from  the  fact  that  for  transmission  scattering,  the  energy  flow 
is  in  the  -z  direction.  This  work  always  assumes  that  the  incident  medium  is  lossless  where 
^{£1+2}  =  0  where  Z{Q}  denotes  the  imaginary  part  of  Q  and  9?{e£+2}  ^  1.  Also,  when 
transmission  scattering  is  considered,  the  substrate  medium  is  likewise  assumed  to  be  lossless 
with  3  {Ej}  =  0  and  ^  1 .  In  general,  however,  Ej  is  complex  with  unrestricted  real  and 

imaginary  parts  such  that  SK{Ei}  can  be  greater  than  or  less  than  zero  and  3{Ei}  is  greater 

than  zero.  Equations  5.3  represent  energy  flow  as  an  integral  sum  over  all  possible  scattered 
wave  vectors  k.  However,  not  all  wave  vectors  contribute  to  scattering,  since  in  Equation 
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5.3a  there  is  no  energy  flow  in  the  z  direction  when  3 {<71+2}  ^  0-  Likewise,  in  Equation 
5.3b,  there  is  no  energy  flow  in  the  — z  direction  when  3{^i'^}  ^  0.  In  considering  reflection 
scattering,  it  follows  that  3{q'ii2}  =  0  unless  k  >  (co/c)'\/ei,+2.  Similarly,  for  transmission 
scattering  when  k  >  we  see  that  3{^f * 0,  These  wave  vector  values  k  =  \k\, 

which  cause  3  {ql +2}  ^  0  and/or  3  P}  *  0,  correspond  to  evanescent  scattered  modes  such 

as  surface  plasmons  or  guided  waves  that  propagate  parallel  to  the  interfaces  and  thus  do 
not  carry  energy  away  from  the  surface.  In  this  work  we  are  interested  only  in  scattered 
fields  that  carry  energy  away  from  the  multilayer  stack,  and  this  can  be  accomplished  in  an 
ad  hoc  fashion  for  Equation  5.3a  (reflection  scattering)  if  we  let 


k^  =  (co/c  )'V  £/, + 2  sin  6  cos  (j) 

and 

(5.4a) 

ky  =  ((o/c)-vei, +2  sin  0  sin  ^ 

(5.4b) 

from  which  it  follows  that  k  =  (co/c)'vei,+2sinO  and 

d^k  -  ((n/c)^e£,+icos0 sin0d0d(|)  =  (a)/c)\+2cos0df2 

(5.4c) 

Similarly  for  Equation  5.3b  (transmission  scattering)  we  let 

k^  =  (oy  c  )'\/ei^  sin  0  cos  <() 
and 

(5.5a) 

ky  =  (co/c  )Vei  sin  0  sin  ([) 

(5.5b) 

from  which  it  follows  that  k  =  (ca/c)Vei^sin0  and 

d^k  —  (Qi/c)^ei  cos  0  sin  QdQd^  =  cos  0dQ 

(5.5c) 

It  is  not  necessary  to  choose  the  and  ky  as  given  in  Equations  5.4  and  5.5  in  an  ad  hoc 

manner.  These  results  also  are  obtained  by  integrating  Equations  4.1  by  the  method  of 
stationary  phase  and  then  proceeding  with  Equation  5.2. 

From  Equation  5.3a  with  Equations  5.4,  we  can  show  that  the  ARS  result  for  scattering 
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into  the  reflection  hemisphere  (medium  £1+2)  into  direction  (0,<t))  is 
1  /dPi^A  (0Vc)Vri2ei42Cose((|flri2l^)+<jgl*l2l^} 

Po\  dCl  /  ^ 

where  the  first  and  second  terms  in  the  {•  •  •}  are  for  p-  and  s-polarization,  respectively,  of 
the  scattered  fieldrelative  to  the  plane  of  scattering  given  by  (k,  z).  The  corresponding  result 
for  scattering  into  the  transmission  hemisphere  (medium  Ej)  is 

1  /^i\  (o^cf<7i\cose((|hn 

PoKdCl/  ^ 

The  angle  brackets  (...)  denote  an  ensemble  average.  The  ensemble  average  is  needed 
because  the  aj^l  2,  and  gf 1 2  depend  on  statistical  properties  of  the  roughness  functions 

Azy(p),  which  are  treated  as  random  variables  for  interfaces  j  =  1  -^L  +  i.  These  statistical 

properties  appear  in  the  power  spectral  density  of  the  interface  roughness,  and  this  is  dis¬ 
cussed  in  the  following  section. 


6.  POWER  SPECTRAL  DENSITY  FUNCTIONS 


Since  the  surface  roughness  function  Azy(p)  for  interface  j  can  never  be  known 

deterministically,  a  statistical  treatment  is  required  to  obtain  an  average  or  expected  result 
for  the  ARS.  In  this  work,  we  assume  that  the  surface  area  A  and  the  statistical  process 

is  stationary,  ergodic,  and  isotropic.  This  means  that  statistical  properties  are  independent 
of  position  and  that  an  ensemble  average  and  spatial  average  yield  identical  results  inde¬ 
pendent  of  direction.  The  averaging  acts  on  the  scattering  coefficients  as  shown  in  Equations 
5.6,  and  these  coefficients  are  rewritten  here  as 


«  =  1  in  =  1 


(6.1a) 


n  =  1  (n  =  1 


(6.1h) 
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,  L  +  lL  +  l  ,  _  _  ,  — 

<  I  >=  S  S  fiX<AZ„(fco  -  -k)) 

Asimsl 

A  L  +  lL  +  l  — •  —4  ^  . 

<l^'!2l  >=  £  2  F.F;<AZ,«:„-i)AZ,(t.-t)> 

l|:=lm  =  l 


(6.1c) 


(6.W) 


where  *  denotes  complex  conjugate.  The  Aj,  Gj,  Bj,  and  Fj  are  given  by 


1  [  (0  c  )^[ej  - 1) _ 1_  v(o)n(»->  -  *)] 


(6.2a) 


(a)/c)Xf^^-'-‘^ 


(62b) 


B,  =  - 


1 

1  (ciVc)V^ 

(6.2c) 


p  _  (i,t + 1)^(1  J  - 1)  _  ^(u + - 1)  J 


(6.2d) 


and  we  define 


-  -  (/SZ,(ko-k)^in(ko-k)) 

g^(K  -k)= - Z - 


(6.3) 


which  is  the  average  cross-power  spectral  density  for  the  spatial  ftequencies  of  the  roughness 
between  interfaces  n  and  m.  This  average  is  easily  shown  to  be  related  to  the  roughness 
correlation  function  by  using  Equation  2.8  in  Equation  6.3,  which  yields 


gnr.(ko-k)  =  j rf'P  j +  P)^«(P) 


(6.4) 


where  the  correlation  function 
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C„(P)  =  <A2,(p  + P)A2,(p)>=  Umj  J  J  d=pAz.(p+ P)Az„(p)|  (6.5) 

gives  the  average  relationship  between  the  roughness  functions  Ar,  and  Az„  for  interfaces 
n  and  m  that  are  laterally  separated  by  lag  length  p.  Note  that  in  Equations  6.4  and  6.5,  we 
have  invoked  the  ergodic,  stationary,  and  isotropic  assumptions  where  G,„„(P)  =  G,„„(P)  and 

p=ia. 


We  consider  two  separate  cases  regarding  the  correlation  of  the  interface  roughness 
between  interfaces  and  these  lead  to  different  ARS  results.  In  the  first  case,  the  roughness 
is  literally  identical  at  each  interface  and  this  is  called  a  correlated  multilayer  stack.  In  the 
second  case,  the  roughness  at  a  given  interface  is  statistically  independent  of  the  roughness 
at  any  other  different  interface  and  this  is  called  an  uncorrelated  multilayer  stack.  The 
statistical  relationship  of  the  roughness  between  interfaces  n  and  m  is  given  by  their 
cross-correlation  function  G*„(P)  where  n^m  or  autocorrelation  function  G„(p)  when 
m—n. 


We  assume  a  functional  form  for  the  correlation  function  G*„(P)  and  calculate  the 

corresponding  power  spectral  density  -  k)  from  Equation  6.4.  In  this  study  we  assume 

for  all  interfaces  n  and  m  that  the  correlation  functions  for  the  correlated  model  can  be 
written  as 

G,„,.(P)  =  G(P)  (6.6a) 

and  this  is  rather  obvious  since  all  interfaces  are  identical.  Thus,  for  this  case,  the  auto-  and 
cross-correlation  functions  are  the  same  at  and  between  all  interfaces.  In  the  case  of  the 
uncorrelated  model,  we  assume 

G»„(P)  =  5„„,G(P)  (6.6b) 

where  the  correlation  function  G(P)  is  the  same  in  both  models.  The  only  difference  is  that 
the  cross-correlation  function  vanishes  in  the  unconelated  model.  The  correlated  and 
uncorrelated  models  are  perhaps  the  simplest  to  consider,  however,  more  complicated 
modeling  can  be  done,  such  as  varying  the  degree  of  correlation  between  interfaces. 

We  write  the  correlation  function  G(P)  as  a  sum  of  an  exponential  and  a  Gaussian  as 
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Gm  =  8le 


-(P/o/ 


(6.7) 


which,  when  used  in  Equation  6.4,  yields  the  power  spectral  density  function  g  (ko  -  k)  as 


g(ko-k)  = 


Inbldi 


[l  +  \ko-kfdl] 


3« 


+  TcSjO^C 


(6.8fl) 


where 


\kQ- k\^  =  k^  +  k^ -2kjc  cos^ 


(6.86) 


The  L  and  S  subscripts  refer  to  long-  and  short-range  roughness  scale,  respectively.  The  5^ 
and  Oi  are  the  long-range  rms  roughness  and  correlation  length,  respectively,  as  the  85  and 
Oj  are  the  short-range  rms  roughness  and  correlation  length,  respectively.  With  these 
assumptions  the  ARS  formulas  for  the  correlated  model  become 


\_  {(i)lcfq2l2^i^2<^osQ  \ 

Po\  dQ  /  (2nM%  1 


L  +  l 

JiAj 


|2 

+ 


L  +  1 

iGj 
>=i  ^ 


(6.9) 


1  /dPf°’^\  (a)/c)^^i%cos0 
da  / 


Z  +  l 

2 

i+1  2 

7  =  1 

+ 

7  =  1 

(2Kfq[% 

The  corresponding  results  for  the  uncorrelated  model  become 

1  \  (oyc)"^ri2ei+2cos  e  ,  ,2. 

p\-dn-/ - '  b?.''''' 


giko-k) 


giko-k) 


(6.10) 


(6.11) 


(co/c)Vf'^eiCos0  fi+i  ,  L+i 

ZI5/+ZIF/ 
U=1  >=1  . 


(2jt)^<7fi2 


gik,-k) 


(6.12) 


where  it  is  emphasized  again  that  for  transmission  scattering  we  assume  a  priori  that 
3{ei}  =  0  and  9i{ei}  >  0. 
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7.  COMPUTER  PROGRAM  AND  NUMERICAL  EXAMPLES 


This  section  describes  how  to  use  a  computer  program  written  to  obtain  numerical 
results  from  the  preceding  ARS  formulas  and  discusses  specific  numerical  examples  of  input 
and  output  data.  Appendix  A  lists  corresponding  symbols  used  in  the  computer  program 
and  symbols  used  in  the  theory  as  given  in  the  preceding  sections.  The  computer  code  is 

given  in  Appendix  B .  In  the  computer  algorithm,  ihe,k^,k,q f\  and  qf^  have  been  normalized 

to  (ale  throughout  This  is  important  to  remember  when  trying  to  compare  the  previously 
outlined  theory  with  the  code.  The  computer  code  consists  of  a  main  program  and  one 
subroutine  called  "ZERO."  The  subroutine  21ERO  computes  the  zero-order  field. 

Taking  into  account  normalization  by  co/c  as  mentioned  above,  the  ARS  for  reflection 
and  transmission  scattering  as  in  Equations  6.9-6.12  are  now  rewritten  as 


1  /dPf+^\  ^t+2£z,+2COS0  f  i  +  t  2  £.  +  ,  2'|  _  _ 


1  ^f^ejcosef  t+i  2  i,+i  2]  _  _ 


1  ^Li2ei.42Cos0f/.+i  , 


VNCORR 


^/ripr 

Po\  da 


?f“e,cos9fi.*i 


where  oVc  =  2n/k  has  been  used.  In  each  of  these  four  equations,  the  first  summation  in  the 
{•  •  •}  brackets  is  the  p-polarized  component  of  the  scattered  field,  whereas  the  second 
summation  is  the  s-polarized  field.  Relating  these  formulas  to  the  corresponding  symbols 
in  the  computer  code  for  p-polarized  reflected  and  transmitted  scattered  light  yields 

1  qij2^L+2COsd  ^-+1  ^ 
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Po\  da  / 


=  RP2(m)  = 


xWlU 


iiA,rg(*o-^) 


(7.6) 


1./E1 

Po 


CORR\  ^}‘^eiCOS0 


da 


)=TPl(m)  =  - 


£+1 


UNCORR' 

1 


_i.M 

p\  da 


\  cos  6^+1  ,  _  — 

)=  TP2  (m)  =  -  Yi  —  £  I  B/s  »o  -  *) 


(7.7) 


(7.8) 


The  corresponding  results  for  s-polarized  reflected  and  transmitted  scattered  light  are 


l/dP“f\_„  «i‘lA.,cose 


£  +  1 


g{K-k) 


(7.9) 


1  /dPy^9°^\  fl'L+sez.+iCosO^+i  ,  _  _ 


(7.10) 


1  /dPy°^  \  cos  0 


£+1 


gih-k) 


(7.11) 


1  ^l'f^eiCos0i+»  2 

)=r52(m)=  ' 


T’oX  da 


mi2  >=> 


I|F/g(^o-^) 


(7.12) 


The  /?Pi I'm),  RP2(m),  etc.,  are  1-D  arrays  for  ARS  at  polar  angle  X(m)  =  THETA  =  0  where 
the  integer  m  =  1  NM  and  NM  is  the  number  of  data  points  calculated.  The  R  and  T 
denote  reflectance  and  transmitted  ARS,  respectively.  The  P  and  S  refer  to  the  scattered 
field  polarization  relative  to  the  plane  of  scattering,  respectively.  Finally,  The  i  and  2  signify 
that  the  multilayer  stack  is  correlated  or  uncorrelated,  respectively. 

Some  of  the  following  statements  reflect  the  fact  that  this  FORTRAN  code  was  run 
on  the  Naval  Air  Warfare  Center  Weapons  Division’s  Cray  supercomputer.  Of  course,  a 
Cray  is  not  needed  since  a  VAX  or  PC  with  the  necessary  memory  will  suffice.  To  run  this 
code  requires  input  of  data,  which  are  read  from  a  file  called  nawc.in  along  with  consistent 
data  in  two  PARAMETER  statements.  The  PARAMETER  statements  are  discussed  in  the 
next  paragraph.  There  are  four  READ  statements.  The  first  statement  reads,  in  degrees,  the 
beginning  (TS),  ending  (TF),  and  increment  (DT)  on  scattering  angle  THETA.  From  these 
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three  values,  the  number  of  data  points  NM  is  calculated.  The  second  statement  reads 
THETAO  (polar  angle  of  incidence  in  degrees),  LAMBDA  (wavelength),  PHIO  (azimuth 
angle  defining  scattering  plane,  in  degrees),  and  L  (number  of  layers  in  the  multilayer  stack). 
The  third  statement  reads  the  rms  roughness  and  correlation  length  values: 
RMSL,SIGL,RMSS,SIGS.  The  last  READ  statement  reads  in  the  multilayer  stack  design, 
layer  by  layer,  as  described  by  their  complex  dielectric  constant  and  physical  thickness. 
These  data  are  read  starting  with  the  substrate,  through  the  L  layers  in  the  stack,  and  ending 
with  the  superstrate.  The  program  prompts  the  user  for  the  unit  number  of  the  output  data 
file.  For  example,  if  you  enter  "35"  as  the  unit  number,  the  output  data  is  stored  in  "fort.35." 

Sample  input  and  output  data  files  are  given  in  Appendixes  C  through  F.  Also  shown 
in  each  appendix  are  sample  plots  of  ARS  versus  polar  scattering  angle  THETA.  The  output 
data  file  lists  the  polarization  angle  of  the  incident  beam  PSI,  the  angles  PHIO  and  THETAO, 
and  the  wavelength  (LAMBDA).  Following  this  are  ARS  output  data  given  in  rows  preceded 
by  a  P  or  S  which  denote  the  polarization  of  the  scattered  field  for  the  data  in  that  row.  In 
rows  preceded  by  a  P  follows  the  polar  scattering  angle  THETA,  correlated  reflected  ARS 
(RPl),  uncorrelated  reflected  ARS  (RP2),  correlated  transmitted  ARS  (TPl),  and  uncor¬ 
related  transmitted  ARS  (TP2).  Similarly,  in  rows  preceded  by  an  S  follows  THETA, 
correlated  reflected  ARS  (RSI),  uncorrelated  reflected  ARS  (RS2),  correlated  transmitted 
ARS  (TSl),  and  uncorrelated  transmitted  ARS  (TS2).  Also  shown  in  each  appendix  is  the 
corresponding  nawc.in  input  data  file.  THETA  values  are  given  as  either  positive  or  negative: 
positive  values  are  for  azimuth  scattering  angle  PHI  =  PHIO;  negative  values  are  for  azimuth 
scattering  angle  PHI  =  PHIO  +180  (degrees). 

All  numerical  values  given  here  used  the  same  rms  and  correlation  length  values  of 
RMSS  =  0.00049  microns,  SIGS  =  0.2  microns,  RMSL  =  0.0005  microns,  and  SIGL  =  2.0 
microns. 

Appendix  C  provides  numerical  input  and  output  results  for  a  28-layer  (L  =  28)  stack, 
angle  of  incidence  THETAO  =  45.0  degrees,  azimuth  angle  PHIO  =  0.0  degrees,  and 
wavelength  LAMBDA  =  0.65  microns.  The  incident  and  substrate  media  dielectric  constants 
are  E(30)  =  (1. 0,0.0)  and  E(l)  =  (2.3104,0.0),  respectively.  The  complete  design  of  the 
multilayer  system  is  given  in  Appendix  C.  Shown  in  Figures  C-l(a)  through  (d)  are  plots 
of  reflected  and  transmitted  ARS  versus  THETA  for  various  polarizations  and  correlated 
or  uncorrelated  interface  roughness. 

Appendix  D  provides  an  example  of  input  and  output  data  for  ARS  from  a  bare  (L  = 
0)  surface  that  includes  total  internal  reflection.  The  incident  and  substrate  media  have 
dielectric  constants  E(2)  =  (4.41,0.0)  and  E(l)  =  (1.0,0.0),  respectively.  The  wavelength 
LAMBDA  =  0.633  microns,  azimuth  angle  PHIO  =  0.0  degrees,  and  angle  of  incidence 
THETAO  =  35.0  degrees,  which  is  beyond  the  critical  angle  of  28.44  degrees.  In  this  case 
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there  is  no  specular  transmitted  beam  but  there  is  transmitted  scattered  light.  Figures  D- 1  (a) 
and  (b)  show  reflected  and  transmitted  ARS  versus  THETA  for  various  polarizations.  Since 
there  is  only  one  interface,  the  correlated  and  uncorrelated  models  give  identical  results. 

Appendix  E  lists  numerical  input  and  output  data  for  an  example  of  out-of-plane  ARS 
where  the  azimuth  angle  PHIO  =  35.0  degrees,  angle  of  incidence  THETAO  =  60.0  degrees, 
and  wavelength  LAMBDA  =  0.633  microns.  The  stack  has  five  layers  (L  =  5)  and  an 
absorbing  metallic  substrate.  In  this  case  there  is  no  transmitted  scattered  light,  and  the 
multilayer  design  is  given  in  the  appendix.  Figures  E-l(a)  through  (d)  plot  reflected  ARS 
versus  THETA  for  all  permutations  of  p-  and  s-polarization  of  the  incident  and  scattered 
light.  Also,  ARS  for  both  correlated  and  uncorrelated  roughness  is  shown. 

Appendix  F  provides  a  final  example  of  input  and  output  data  for  a  single  layer  case 
(L  =  1)  where  the  incident  medium  is  E(3)  =  (4.41,0.0),  followed  by  a  partially  opaque 
0.025-micron-thick  metal  film  with  dielectric  constant  E(2)  =  (-16.0,0.5)  and  an  air  substrate 
with  E(  1 )  =  ( 1 .0,0.0).  The  wavelength  LAMBDA =0.633  microns,  THETAO =29.6  degrees, 
and  PHIO  =  0.0  degrees.  This  case  is  of  interest  because,  for  p-polarized  incident  light  at 
this  angle  of  incidence,  direct  coupling  to  a  surface  plasma  resonance  in  the  metal  film  does 
occur  and  this  has  an  effect  on  the  resulting  scattered  light.  Figures  F-l(a)  through  (d)  show 
plots  of  reflected  and  transmitted  ARS  versus  THETA  for  p-  and  s-polarization  of  the  incident 
and  scattered  fields. 
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Appendix  A 

CORRESPONDENCE  BETWEEN  SYMBOLS  USED  IN  THE 
COMPUTER  CODE  AND  THEORY 


Listed  below  are  the  symbol  relationships  between  the  computer  listing  in  Appendix  B 
and  the  theory  as  discussed  in  the  main  body  of  this  publication. 
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Appendix  B 

FORTRAN  LISTING  OF  MULTILAYER  SCATTERING  SOURCE  CODE 


C  TfflS  PROGRAM  IS  DESIGNED  TO  PREDICT  ANGLE  RESOLVED 
C  SCATTERING  (ARS)  DUE  TO  ROUGHNESS  PRESENT  AT  THE  INTERFACES 
C  OF  MULTILAYERED  OPTICAL  SURFACES.  THE  SCATTERING  IS 
C  CALCULATED  AS  A  REFLECTED  OR  TRANSMITTED  FIELD,  ms 
C  VALID  ONLY  WHEN  THE  RMS  ROUGHNESS  VALUES  OF  THE 
C  INTERFACES  ARE  MUCH  LESS  THAN  THE  INCIDENT  WAVELENGTH. 

C  THE  PROGRAM  AS  WRITTEN  HERE  CONSIDERS  THE  CASES  OF  (1) 

C  IDENTICAL  ROUGHNESS  AT  EACH  INTERFACE  (CORRELATED), 

C  (2)  INDEPENDENT  ROUGHNESS  BETWEEN  DIFFERENT  INTERFACES, 

C  (UNCORRELATED),  AND  (3)  A  PARTIALLY  CORRELATED  MODEL.  THE 
C  PARTIAL  CORRELATION  MODEL  ASSUMES  THAT  THE  LONG-RANGE- 
C  SCALE  ROUGHNESS  IS  CORRELATED  BETWEEN  INTERFACES  AND 
C  THE  SHORT-RANGE-SCALE  ROUGHNESS  IS  UNCORRELATED  BETWEEN 
C  INTERFACES.  THE  OUTPUT  IS  THE  ARS,  i.e.,  SCATTERED  POWER 
C  PER  UNIT  SOLID  ANGLE  NORMALIZED  TO  THE  INCIDENT  POWER 
C  VERSUS  POLAR  ANGLE  OF  SCATTERING.  FOR  REFLECTION 
C  SCATTERING,  THE  POLAR  SCATTERING  ANGLE  IS  IN  THE  INCIDENT 
C  MEDIUM.  FOR  TRANSMISSION  SCATTERING.  THE  POLAR  SCATTERING 
C  ANGLE  IS  IN  THE  SUBSTRATE  MEDIUM. 

C 

C  NL  =  L  +  2  (SEE  MAIN  PROGRAM  AND  SUBROUTINE  "ZERO") 

C  NP  =  MAXIMUM  NUMBER  OF  OUTPUT  DATA  POINTS  TO  CALCULATE 
C 

C  IN  THE  PARAMETER  STATEMENT  OF  THE  FOLLOWING  PROGRAM  LISTING, 
C  THE  NUMBER  OF  LAYERS,  L  =  5,  AND  THUS  NL  =  7. 


PARAMETER(NP  =  400 J^L  =  7) 

DIMENSION  X(NP) 

COMPLEX  XI1(0:NL) ,  XI2(0:NL) ,  MU1(0:NL) .  MU2(0:NL) 
COMPLEX  ETA1(0:NL) ,  ETA2(0:NL) 

COMPLEX  ETA3(0:NL) ,  ETA4(0:NL) 

COMPLEX  PI  1(0:NL)  ,  P12(0:NL)  ,  P21(0:NL) ,  P22(0:NL) 
COMPLEX  S11(0;NL) ,  S12(0:NL) ,  S21(0:NL) ,  S22(0:NL) 
COMPLEX  E(NL) ,  Q0(NL)  ,  Q1(NL) 

COMPLEX  PP ,  PM  ,  SP ,  SM  ,  Z1 1 .  Z12  ,  Z21 ,  Z22 
COMPLEXI,Z,Zl  ,Z2,Z3,Z4,ZU.ZB,ZE 
COMPLEX  ZUMl ,  ZUM2 ,  ZUM3  .  ZUM4 
COMPLEX  AJ  ,  BJ  ,  FJ  ,  GJ 
COMPLEX  Yll ,  Y12  ,  Y21 ,  Y22 
COMPLEX  XI 1 ,  X12  ,  X21 ,  X22 
COMPLEX  Y1(NL)  ,  Y2(NL)  ,  X0(NL) 

COMPLEX  E0X(NL)  ,  E0Z(NL)  ,  E0Y(NL) 

REAL  LAMBDA,  K,K0 
REALD(NL),T(NL) 

REAL  SUMl ,  SUM2 

REAL  RP1(NP) ,  RS1(NP)  ,  RP2(NP) 

REAL  RS2(NP) ,  RP3(NP) ,  RS3(NP) 
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REAL  TP1(NP)  ,  TS1(NP) ,  TP2(NP) 

REAL  TS2(NP) ,  TP3(NP) .  TS3(NP) 

INTEGER  NM 

OPEN(unit  =  20,  file  =  ’nawc.in’ ,  status  =  ’old’) 

READ(20.*)  TS,TF,DT 
READ(20,*)  THETAO  J.AMBDA  PfflOX 
READ(20,*)  RMSL,SIGLJIMSS.SIGS 

PI  =  ACOS(-l.O) 

L1=L+1 
L2  =  L  +  2 

WC  =  2*PI/LAMBDA 
I  =  (0.0, 1.0) 

DR  =  PI/180. 

READ(20,*)  (E(M),T(M),  M  =  1,  L2) 

C  INPUT  PARAMETERS: 

C  TS  =  START  VALUE  (DEG)  OF  POLAR  SCATTERING  ANGLE  THETA 
C  TF  =  FINAL  VALUE  (DEG)  OF  POLAR  SCATTERING  ANGLE  THETA 
C  DT  =  INCREMENT  ON  POLAR  ANGLE  THETA  (DEGREES) 

C  THETAO  =  ANGLE  OF  INQDENCE  (DEG)  MEASURED  FROM  THE  NORMAL 
C  PmO  =  AZIMUTHAL  SCATTERING  ANGLE  DESIRED  (DEG).  FOR 
C  EXAMPLE,  IF  YOU  WANT  TO  LOOK  AT  SCATTERING  IN  THE  PLANE  OF 
C  INCIDENCE,  SET  PHIO  =  0.  THIS  PROGRAM  AUTOMATICALLY  SHIFTS 
C  PHIO  180  DEGREES  FOR  SCATTERING  IN  THE  QUADRANT  CONTAINING 
C  THE  INCOMING  BEAM. 

C  LAMBDA  =  INCIDENT  WAVELENGTH 

C  L  =  INTEGER  NUMBER  OF  LAYERS  IN  MULTILAYER  STACK 

C 

C  E(J)  =  COMPLEX  DIELECTRIC  CONSTANT  ARRAY 
C 

C  E(l)  =  SUBSTRATE 
C  E(2)  =  FIRST  LAYER 
C  . 

C  . 

C  . 

C  E(L+1)  =  OUTERMOST  LAYER 

C  E(L+2)  =  REGION  ABOVE  STACK  -  SUPERSTRATE 

C 

C  T(J):  THICKNESSES  OF  THE  ITH  LAYER 
C 

C  T(l)  =  NOT  USED 

C  T(2)  PHYSICAL  THICKNESS  OF  LAYER  NEXT  TO  SUBSTRATE 
C  . 

C  . 

C  . 

C  T(L+1)  =  OUTERMOST  LAYER  PHYSICAL  THICKNESS 
C  T(L+2)  =  NOT  USED 
C 

D(1)  =  0.0 

C  D(J)  =  ABSOLUTE  THICKNESS  OF  MULTILAYER  STACK  AT  INTERFACE  J 

IF(L.GT.0)  THEN 
DO  706  J  =  2  ,  LI 
706  D(J)  =  D(J-1)  +  T(J) 

ELSE  IF(L.EQ.0)  THEN 
ENDIF 

C  NM.LEJNIP  IS  THE  CALCULATED  NUMBER  OF  OUTPUT  DATA  POINTS 
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XM  =  (TF-TS)/DT+L0 
NM  =  INT(XM)+  1 

1001  CONTINUE 

KO  =  SQRT(REAL(E(L2)))*SIN(DR*THETA0) 

DO301J=l.L2 

301  Q0(J)  =  CSQRT(E(J)  -  K0**2) 

CALL  ZERO(EOX^OZ£OY.QO£L)  JC0JL.WC) 

NU  =  2 

IF(REAL(E(1)).LE.0.0)  NU  =  1 
D0121JJ=1.2 

C  TWO  VALUES  OF  POLARIZATION  ANGLE  SIGMA  ARE  USED,  BUT 
C  ANY  RANGE  OF  VALUES  (IN  RADIANS)  MAY  BE  USED 
C 

C  WHEN  JJ  =  1  ->  P-POLARIZED  INCIDENT  BEAM 
C  WHEN  JJ  =  2  ->  S-POLARIZED  INCIDENT  BEAM 
C 

SIGMA  =  (JJ-l)*(PI/2.) 

POC  =  COS(SIGMA) 

POS  =  SIN(SIGMA) 

D0122U=1.ND 

C  WHEN  D  =  1  ->  REFLECTION  SCATTERING 
C  WHEN  U  =  2  ♦>  TRANSMISSION  SCATTERING 
C  IF(REAL(E(1))LE.0.0)  TRANSMISSION  SCATTERING  IS  OMITTED 
C 

DO120n=l,NM 

C  THETA  ==  POLAR  ANGLE  OF  SCATTERING 

THETA  =TS  +  (n-  1)*DT 
X(II)  =  THETA 

C  NOTE  DIFFERENT  K  FOR  REFLECTION  (D=l)  AND  TRANSMISSION  (U=2)  SCATTERING 
IF(U.EQ.l)  THEN 

K  =  SQRT(REAL(E(L2)))*ABS(SrN(DR*THETA)) 

ELSE  IF(U.EQ.2)  THEN 

K  =  SQRT(REAL(E(1)))*ABS(SIN(DR*THETA)) 

ENDIF 

C  PHI  =  AZIMUTHAL  ANGLE  OF  SCATTERING 
C  IF  SIGN  OF  THETA  POSITIVE,  PHI  =  PHIO 
C  IF  SIGN  OF  THETA  NEGATIVE,  PHI  =  PHIO  +  180 

IFCTHETALTO.)  PHI  =  DR*(PHI0  +  180.) 

IF(THETA.GE.O.)  PHI  =  DR*PHI0 

AZC  =  COS(PHI) 

AZS  =  SIN(PHI) 

DO  310  J=  1  ,L2 

310  Q1(J)  =  CSQRT(E(J)  -  K**2) 

C  SOURCE  TERMS 
DO302J=l,Ll 
ZE=:E(J+1)-E(J) 

Y1(J)  =  I*WC*ZE*K*E0Z(J)*POC/E(J+l) 

Y2(J)  =  I*WC*ZE*(E0X(J)*AZC*POC  +  E0Y(J)*AZS*POS) 

X0(J)  =  I*WC*ZE*(E0X(J)*AZS*POC  -  E0Y(J)*AZC*POS) 
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302  CONTINUE 

Yll  =  (1..0.) 

Y12  =  (0..0.) 

Y21  =  (0.,0.) 

Y22  =  (l.,0.) 

X11  =  (1.,0.) 

X12  =  (0.,0.) 

X21  =  (0..0.) 

X22  =  (l.,0.) 

D0311M=  1  .LI 

ZB  =  I*WC*(Q1(M+1)  -  Q1(M))*D(M) 

ZU  =  I*WC*(Q1(M+1)  +  Q1(M))*D(M) 

Z  =  CSQRT(E(M)*E(M+1)) 

PP  =  (Q1(M+1)*E(M)  +  Q1(M)*E(M+1))/(2*Q1(M)*Z) 
PM  =  (Q1(M+1)*E(M)  -  Q1(M)*E(M+1))/(2*Q1(M)*Z) 

Z11  =  PP*CEXP(ZB) 

Z12  =  PM*CEXP(-ZU) 

Z21=PM*CEXP(ZU) 

Z22  =  PP*CEXP(-ZB) 

P11(M)  =  Y11*Z11  +  Y12*Z21 
P12(M)  =  Y1 1*Z12  +  Y12*Z22 
P21(M)  =  Y21*Z11  +  Y22*Z21 
P22(M)  =  Y21*Z12  +  Y22*Z22 

Y11  =  P11(M) 

Y12  =  P12(M) 

Y21=P21(M) 

Y22  =  P22(M) 

SP  =  (Q1(M)  +  Q1(M+1))/(2*Q1(M)) 

SM  =  (Q1(M)  -  Q1(M+1))/(2*Q1(M)) 

Zll  =  SP*CEXP(ZB) 

Z12  =  SM*CEXP(-ZIJ) 

Z21  =  SM*CEXP(ZU) 

Z22  =  SP*CEXP(-ZB) 

S11(M)  =  X11*Z11  +  X12*Z21 
S12(M)  =  X11*Z12  +  X12*Z22 
S21(M)  =  X21*Z11  +  X22*Z21 
S22(M)  =  X21*Z12  +  X22*Z22 

X11  =  S11(M) 

X12  =  S12(M) 

X21  =  S21(M) 

X22  =  S22(M) 

311  CONTINUE 

P11(0)  =  (1.,0.) 

P12(0)  =  (0.,0.) 

P21(0)  =  (0.,0.) 

P22(0)  =  (l.,0.) 

S11(0)  =  (1.,0.) 

S12(0)  =  (0.,0.) 

S21(0)  =  (0.,0.) 

S22(0)  =  (l.,0.) 

DO  3001  =  0, L 
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Z1  =  CEXP(-I*WC*Ql(J+irD(J+l)) 

Z2  =  CEXP(  I*WC*Q1(J+1)*D(J+1)) 

ETA1(J)  =  P11(J)*Z1  +  P12(J)*Z2 
ETA2(J)  =  P11(J)*Z1  -  P12(J)*Z2 
ETA3(J)  =  P21(J)*Z1  +  P22(J)*Z2 
ETA4(J)  =  P21(J)*Z1  -  P22(J)*Z2 

XI1(J)  =  S11(J)*Z1  -  S12(J)*Z2 
Xn(J)  =  S21(J)*Z1  -  S22(J)*Z2 

MU1(J)  =  P11(L1)*ETA3(J)  -  P21(L1)*ETA1(J) 

MU2(J)  =  P11(L1)*ETA4(J)  -  P21(L1)*ETA2(J) 

300  CONTINUE 

ARC  =  K**2  +  K0**2  -  2*K*K0*AZC 
IF(ARGJLE.O.)  DELK  =  0. 

IF(ARG.GT.O.)  DELK  =  SQRT(ARG) 

DELK  =  WC*DELK 
C 

C  THIS  PROGRAM  ASSUMES  A  TWO  PART  AUTOCORRELATION  FUNCTION 
C  (ACF): 

C  ACF  =  EXPONENTIAL  +  GAUSSIAN 
C 

C  THE  RESULTING  SPECTRAL  DENSITY  FUNCTION  NEEDS  FOUR  (4)  INPUT 
C  PARAMETERS. 

C 

C  RMSL  =  LONG-RANGE  ROUGHNESS  RMS  VALUE 
C  RMSS  =  SHORT-RANGE  ROUGHNESS  RMS  VALUE 
C  SIGL  =  LONG-RANGE  ROUGHNESS  CORRELATION  LENGTH 
C  SIGS  =  SHORT-RANGE  ROUGHNESS  CORRELATION  LENGTH 
C 

GKL  =  2*PI*(RMSL*SIGL)**2/((1.+(SIGL*DELK)**2)**1.5) 

GKS  =  PI*(RMSS*SIGS)**2*EXP(-(SIGS*DELK/2.)**2) 

GK  =  GKS  +  GKL 

ZUMl  =  (O.O.O.O) 

ZUM2=  (0.0,0.0) 

ZUM3  =  (0.0, 0.0) 

ZUM4=  (0.0,0.0) 

SUM1=  0.0 
SUM2=  0.0 
SUM3=  0.0 
SUM4=  0.0 

D0313J=1,L1 

C  ZUMl  =  SUM  FOR  P-POL  SCATTERED.  CORRELATED,  REFLECTED 
C  SUMl  =  SUM  FOR  P-POL  SCATTERED,  UNCORRELATED,  REFLECTED 

Z  =  CSQRT(E(J)) 

AJ  =  (Z*Y1(J)*ETA1(J-1)/Q1(J)  -  Y2(J)*ETA2(J-1)/Z)/(2*P11(L1)) 

ZUM1=ZUM1  +  AJ 

SUMl  =  SUMl  -I-  (CABS(AJ))**2 

C  ZUM2  =  SUM  FOR  S-POL  SCATTERED.  CORRELATED.  REFLECTED 
C  SUM2  =  SUM  FOR  S-POL  SCATTERED.  UNCORRELATED,  REFLECTED 

GJ  =  -X0(J)*XI1(J-1)/(2*Q1(J)*S11(L1)) 

ZUM2  =  ZUM2-i-GJ 

SUM2  =  SUM2  +  (CABS(GJ))**2 

C  ZUM3  =  SUM  FOR  P-POL  SCATTERED,  CORRELATED,  TRANSMITTED 
C  SUM3  =  SUM  FOR  P-POL  SCATTERED,  UNCORRELATED,  TRANSMITTED 
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BJ  =  -(Z*Y1(J)*MU1(J-1)/Q1(J)  -  Y2(JrMU2(J-l)/Z)/(2*Pll(Ll)) 

ZUM3  =  ZUM3  +  BJ 

SUM3  =  SUM3  +  (CABS(BJ))**2 

C  ZUM4  =  SUM  FOR  S-POL  SCATTERED.  CORRELATED,  TRANSMITTED 
C  SUM4  =  SUM  FOR  S-POL  SCATTERED,  UNCORRELATED,  TRANSMITTED 

FJ  =  X0{J)*(S11(L1)*XI2(J-1)  -  S21(L1)*XI1(J-1))/(2*Q1(J)*S11(L1)) 

ZUM4  =  ZUM4-hFJ 

SUM4  =  SUM4  +  (CABS(FJ))**2 


313  CONTINUE 

Z  =  (Q1(L2)/QO(L2))*(E(L2)/LAMBDA**2)*COS(DR*THETA) 

CR  =  REAL® 

Z  =  (Q1(1)/QO(L2))*(E(1)/LAMBDA**2)*COS{DR*THETA) 

CT  =  REAL® 

IF(U£Q.l)  THEN 

RPiai)  =  CR*(CABS(ZUM1))**2 
RP2ai)  =  CR*SUM1 

C  RP3  IS  PARTIALLY  CORRELATED  MODEL 
RP3ai)  =  GKL*RPiaD  +  GKS*RP2(II) 

RPiaO  =  GK*RP1(II) 

RP2ai)  =  GK*RP2(II) 

RSiai)  =  CR*(CABS(ZUM2))**2 
RS2ai)  =  CR*SUM2 

C  RS3  IS  PARTIALLY  CORRELATED  MODEL 
RS3(II)  =  GKL*RSiai)  +  GKS*RS2(II) 

RSiai)  =  GK*RSiai) 

RS2(n)  =  GK*RS2ai) 

ELSE  IF(UEQ.2)  THEN 

TP1(II)  =  CT*(CABS(ZUM3))**2 
TP2ai)  =  CT*SUM3 

C  TP3  IS  PARTIALLY  CORRELATED  MODEL 
TP3ai)  =  GKL*TPiai)  +  GKS*TP2(II) 

TPiai)  =  GK*TPiai) 

TP2ai)  =  GK*TP2ai) 

TS1(II)  =  ABS(CT*(CABS(ZUM4))**2) 

TS2ai)  =  ABS(CT*SUM4) 

C  TS3  IS  PARTIALLY  CORRELATED  MODEL 
TS3ai)  =  GKL*TSiaD  +  GKS*TS2ai) 

TSl(II)  =  GK*TSiai) 

TS2ai)  =  GK*TS2(II) 

ENDIF 

120  CONTINUE 
122  CONTINUE 
WRnE(IUNIT,54) 

PRINT*,’  POLARIZATION  ANGLE  SIGMA  =  ’.SIGMA*  180/PI,’  DEGREES’ 
PRINT*,’  DATA  OUTPUT  UNIT  NUMBER  ?  ’ 

READ*TUNIT 

WRnE(IUNIT,55) 

WRITE(IUNIT.50)  SIGMA*  180/PI.LAMBDA 
WRITE(IUNIT,51)  THETAOJPHIO 
WRITE(IUNIT.54) 
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WRITE(IUNIT,11) 

WRITE(IUNIT,12) 

WRITECIUNIT^) 

D0446J  =  1,NM 

WRrrE(IUNIT,9)  X(J)JIP1(J)JIP2(J),TP1(J).TP2(J) 
WRrrE(IUNIT,10)X(J)JlSl(J)JlS2(J),TSl(J),TS2(J) 

446  CONTINUE 

121  CONTINUE 

9  FORMATC  P’ JJ10.2.4E12.4) 

10  FORMATC  S’J?10.2,4E12.4) 

11  FORMATC  MX,‘THETA\5X,’RPr,9X.’RP2’,9X,’TPl’,9X,’TP2’) 

12  FORMATC  MX,’THETA’,5X,’RSr,9X,’RS2’.9X.’TSr,9X,’TS2’) 

50  FORMATC  SIGMA  =’J?10.2,’  LAMBDA  =’J8.4) 

51  FORMATC  THETA0=’JF10.2,’  PfflO  =’f8.2) 

54  FORMATC  ’) 

FORMATC  '*'’*"*"*"*"*“*'’*"i"*‘’*"*‘’*"*‘'*"*"*"*"*"*"*"*‘’*'’*“*‘’*“*"*"*‘’*"*‘’*"*"*"*"*“*"*"*"*‘'*‘’*"*‘’*"*"*‘'*"*'*'*“*‘'*"*‘*’*'’) 
STOP 
END 

SUBROUTINE  ZERO(EOX30ZJSOY,QOJE.DJCOX.WC) 

C  NL  =  NUMBER  OF  LAYERS(L)  +  2  SHOULD  MATCH  MAIN  PROGRAM 
PARAMETER(NL  =  7) 

COMPLEX  QO(NL)X(NL)X0X(NL)X0Z(NL)JEOY(NL) 

COMPLEX  A0(NL)30(NL),G0(NL)JF0(NL) 

COMPLEX  X11,X12,X21,X22 
COMPLEX  Y11,Y12,Y21,Y22 
COMPLEX  S011(0:NL),S012(0:NL) 

COMPLEX  S021(0:NL).S022(0:NL) 

COMPLEX  P011(0:NL)F012(0M.) 

COMPLEX  P021(0:NL)F022(0:NL) 

COMPLEX  PPXM,SP,SMX11X12Z21Z22 
COMPLEX  I2X1X22UXB 
REAL  WC,KOXP,TPXS,TS,D(NL) 

INTEGER  LX1X2 

I  =  (0.0.1.0) 

L1=L+1 
L2  =  L  +  2 


Y11  =  (1.,0.) 
Y12  =  (0.,0.) 
Y21  =  (0.,0.) 
Y22  =  (1..0.) 

Xll  =  (1..0.) 
X12  =  (0.,0.) 
X21  =  (0.,0.) 
X22  =  (1..0.) 


D0311M=1,L1 

ZB  =  I*WC*(Q0(M+1)  -  Q0(M))*D(M) 

ZU  =  PWC*(Q0(M+1)  +  Q0(M))*D(M) 

Z  =  CSQRT(E(M)*E(M+1)) 

PP  =  (Q0(M+1)*E(M)  +  Q0(M)*E(M+1))/(2*Q0(M)*Z) 
PM  =  (Q0(M+1)*E(M)  -  Q0(M)*E(M+1))/(2*Q0(M)*Z) 
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Z11  =  PP*CEXP(ZB) 

Z12  =  PM*CEXP(-ZU) 

Z21  =  PM*CEXP(ZU) 

Z22  =  PP*CEXP(-ZB) 

P011(M)  =  Y11*Z11  +  Y12*Z21 
P012(M)  =  Y11*Z12  +  Y12*Z22 
P021(l^  =  Y21*Z1 1  +  Y22*Z21 
P022(M)  =  Y21*Z12  +  Y22*Z22 

Y11  =  P011(M) 

Y12  =  P012(M) 

Y21  =  P021(M) 

Y22  =  P022(M) 

SP  =  (Q0(M)  +  Q0(M+1))/(2*Q0(M)) 

SM  =  (Q0(M)  -  Q0(M+1))/(2*Q0(M)) 

Zll  =  SP*CEXP(ZB) 

Z12  =  SM*CEXP(-ZU) 

Z21  =  SM*CEXP(ZU) 

Z22  =  SP*CEXP(-ZB) 

S011(M)  =  X11*Z11  +  X12*Z21 
S012(M)  =  X11*Z12  +  X12*Z22 
S021(M)  =  X21*Z1 1  +  X22*Z21 
S022(l^  =  X21*Z12  +  X22*Z22 

X11  =  S011(M) 

X12  =  S012{M) 

X21  =  S021(M) 

X22  =  S022(M) 

311  CONTINUE 

P011(0)  =  (l.,0.) 

P012(0)  =  (0.,0.) 

P021(0)  =  (0..0.) 

P022(0)  =  (1..0.) 

S011(0)  =  (l.,0.) 

S012(0)  =  (0.,0.) 

S021(0)  =  (0.,0.) 

S022(0)  =  (l.,0.) 

D0104J  =  0,L1 

A0(J+1)  =  -Q0(L2)*P012(J)/(P01 1(L1)*Q0(J+1)) 

B0(J+1)  =  Q0{L2)*P011(J)/(P011(L1)*Q0(J+1)) 

F0(J+1)  =  Q0(L2)*S011(J)/(S0ll(LirQ0(J+l)) 

G0(J+1)  =  -Q0(L2)*S012(J)/(S011(L1)*Q0(J+1)) 

104  CONTINUE 

C  INCIDENT  FIELD  COMPONENTS  EVALUATED  AT  D(J) 

DO110J  =  l,Ll 

Z1  =  CEXPa*WC*Q0(J)*D(J)) 

Z2  =  CEXP(-I*WC*Q0{J)*D(J)) 

Z  =  CSQRT(E(J)) 

E0X(J)  =  Q0(J)*(A0(J)*Z1  +  B0(J)*Z2)/Z 
E0Y(J)  =  -G0{J)*Z1  -  F0(J)*Z2 
E0Z(J)  =  -K0*(A0(:)*Z1  -  B0(J)*Z2)yZ 

110  CONTINUE 
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C  SPECULAR  REFLECTANCE 
RP=  (CABS(A0(L2)))**2 
RS  =  (CABS(G0(L2)))**2 
C  SPECULAR  TRANSMISSION 

TP  *  REAL(Q0(1)/Q0(L2))*(CABS(B0(1)))**2 
TS  =  REAL(Q0(1)/Q0(L2))*(CABS(F0(1)))**2 

SUMS  =RS  +  TS 
SUMP  =RP  +  TP 

IF(SUMS.GT.1.001.0R.SUMP.GT.1.001)  THEN 

PRINT*,’  SUM  OF  ZERO-ORDER  ENERGY  GREATER  THAN  UNITY  ’ 

PRINT* ,’SUMP=  ’.SUMP,’  SUMS  =  ’.SUMS 

ELSE 

RETURN 

ENDIF 

RETURN 

END 
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Appendix  C 

INPUT  DATA  AND  OUTPUT  RESULTS  FOR  REFLECTION  AND 
TRANSMISSION  ARS  FROM  A  MULTILAYER  STACK  WITH  28  LAYERS 


OUTPUT  DATA  FE£: 


SIGMA  =  00.00  LAMBDA  =  0.65 

THETA0  =  45.00  PHIO  =  0.00 


THETA 

RPl 

RP2 

TPl 

TSl 

THETA 

RSI 

RS2 

TP2 

TS2 

P  =»  p-POL  scatt.  S  =»  s-POL  scatt. 

P 

-80.00 

0.5510E-06 

0.2279E-06 

0.1392E-06 

0.4680E-06 

S 

-80.00 

0.7422E-34 

0.1363E-33 

0.6056E-35 

0.2224E-34 

INPUT  DATA  FILE: 

P 

-60.00 

0.5691E-05 

0.1488E-05 

0.2453E-07 

0.1772E-05 

s 

-60.00 

0.1993E-33 

0.5344E-33 

0.9153E-33 

0.3348E-33 

-80.,80.,20. 

p 

-40.00 

0.1307E-04 

0.5629E-05 

0.3252E-05 

0.1223E-05 

45.,0.65,0.03 

s 

-40.00 

0.4719E-33 

0.7617E-33 

0.1669E-34 

0.9735E-35 

.0005^.0..00049,0,2 

p 

-20.00 

0.4870E-04 

0.1396E-04 

0.4698E-06 

0.3695E-06 

(2.3104,0.0),0.0 

s 

-20.00 

0,3651E-32 

0.1595E-32 

0.2313E-34 

0.2199E-34 

(4.41.0.0), 0.07738 

p 

0.00 

0.4693E-04 

0.3396E-04 

0.8650E-04 

0.7576E-04 

(2.1316, 0.0), 0.1113 

s 

0.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

(4.41,0.0).0.07738 

p 

20.00 

0.1020E-03 

0.5574E-04 

0.5603E-05 

0.5919E-05 

(2.13 16,0.0).0.1 113 

s 

20.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

(4.41,0.0),0.07738 

p 

40.00 

0.6168E-03 

0.1823E-02 

0.4575E-04 

0.3884E-04 

(2.1316,0.0),0.1113 

s 

40.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

(4.41,0.0),0.30952 

p 

60.00 

0.101^-03 

0.9758E-04 

0.1352E-05 

0.4601E-04 

(2.1316.0.0),0.1113 

s 

60.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

(4.41,0.0).0.07738 

(2.1316.0.0),0.1113 

SIGMA 

=  90.00 

LAMBDA 

0.65 

(4.41.0.0).0.07738 

THETAO 

=  45.00 

PHIO 

=  0.00 

(2.1316,0.0),0.1113 

THETA 

RPl 

RP2 

TPl 

TSl 

(4.41,0.0),0.07738 

THETA 

RSI 

RS2 

TP2 

TS2 

(2.1316,0.0),0.1113 

p 

-80.00 

0.8865E-34 

0.5082E-34 

0.798  lE-35 

0.7769E-34 

(4.41,0.0),0.07738 

s 

-80.00 

0.2955E-06 

0.1567E-05 

0.733  lE-07 

0.1581E-06 

(2.13 16.0.0).0.1 113 

p 

-60.00 

0.6772E-33 

0.3114E-33 

0.4268E-36 

0.3414E-33 

(4.41,0.0),0.07738 

s 

-60.00 

0.7810E-06 

0.6234E-05 

0.6857E-05 

0.2264E-05 

(2.1316.0.0),0.1113 

p 

-40.00 

0.1651E-32 

0.1733E-32 

0.5596E-33 

0.2206E-33 

(4.41.0.0).0.07738 

s 

-40.00 

0.2238E-05 

0.854  lE-05 

0.2812E-07 

0.4702E-07 

(2.1316.0.0),0.1113 

p 

-20.00 

0.7515E-32 

0.3610E-32 

0.4698E-34 

0.5450E-34 

(4.41,0.0).0.30952 

s 

-20.00 

0.1810E-04 

0.1409E-04 

0.3415E-07 

0.5619E-07 

(2.1316,0.0),0.1113 

p 

0.00 

0.1364E-32 

0.9866E-33 

0.2513E-32 

0.2201E-32 

(4.41,0.0),0.07738 

s 

0.00 

0.2432E-04 

0.3003E-04 

0.5186E-04 

0.6718E-04 

(2.1316,0.0),0.1113 

p 

20.00 

0.2965E-32 

0.1619E-32 

0.1628E-33 

0.1720E-33 

(4.41.0.0),0.07738 

s 

20.00 

0.6900E-04 

0.5369E-04 

0.6423E-06 

0.1057E-05 

(2.1316,0.0),0.1113 

p 

40.00 

0.1792E-31 

0.5298E-31 

0.1329E-32 

0.1128E-32 

(4.41,0.0),0.10369 

s 

40.00 

0.6030E-03 

0.2302E-02 

0.8518E-06 

0.1425E-05 

(2.1316,0.0),0.14469 

p 

60.00 

0.2953E-32 

0.2835E-32 

0.3929E-34 

0.1337E-32 

(1.00,0.),0.0 

s 

60.00 

0.5248E-04 
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e.  DEGREES 

(a) 


e,  DEGREES 


(b) 


FIGURE  C-1.  Log[ARS]  Versus  Polar  Scattering  Angle  6.  Plotted  is  scattering  in  the  plane  of  incidence 
from  a  28-layer  multilayer  stack  with  the  incident  beam  of  wavelength  0.65  microns  incident  at  45  degrees. 
Positive  and  negative  0  values  are  for  scattering  azimuth  angles  <(>  of  0  and  180  degrees,  respectively.  The 
solid  and  dashed  lines  refer  to  ARS  when  the  multilayer  stack  interface  roughness  is  correlated  and  uncorrelated, 
respectively.  Reflection  scattering  is  shown  in  (a)  and  (b),  and  transmission  scattering  is  shown  in  (c)  and 
(d).  The  incident  and  scattered  polarization  for  these  four  plots  are  (a)  p  -POL  ^  p  -POL,  (b) 
s -POL  s -POL,  (c)  p -POL  -» p -POL,  and  (d)  s -POL  s -POL. 


40 


NAWCWPNS  TP  8084 


Appendix  D 

INPUT  DATA  AND  OUTPUT  RESULTS  FOR  REFLECTION  AND 
TRANSMISSION  ARS  FROM  AN  UNCOATED  SURFACE  AND 
TOTAL  INTERNAL  REFLECTION  OF  THE  INCIDENT  BEAM 


OUTPUT  DATA  FILE: 

SIGMA  =  00.00  LAMBDA  =  0.633 

THETA0  =  35.00  PHIO  =  0.00 


P  =»  p-POL  scatt.  S  =>  s-POL  scatt. 


THETA 

RPl 

RP2 

TPl 

TSl 

THETA 

RSI 

RS2 

TP2 

TS2 

P 

-80.00 

0.1431E-08 

0.143  lE-08 

0.3313E-07 

0.3313E-07 

S 

-80.00 

0.7522E-36 

0.7522E-36 

0.9796E-36 

0.9796E-36 

P 

-60.00 

0.1473E-07 

0.1473E-07 

0.1589E-06 

0.1589E-06 

INPUT  DATA  FILE; 

S 

-60.00 

0.8622E-35 

0.8622E-35 

0.8856E-35 

0.8856E-35 

P 

-40.00 

0.5331E-07 

0.5331E-07 

0.4229E-06 

0.4229E-06 

-80..80.,20. 

S 

-40.00 

0.5985E-34 

0.5985E-34 

0.3490E-34 

0.3490E-34 

35.0,0.633,0.0,0 

P 

-20.00 

0.2433E-05 

0.2433E-05 

0.1095E-05 

0.1095E-05 

.0005,2.0„00049,0.2 

S 

-20.00 

0.4018E-33 

0.4018E-33 

0.1161E-33 

0.1161E-33 

(1.0,0.0),0.0 

P 

0.00 

0.2470E-04 

0.2470E-04 

0.2667E-05 

0.2667E-05 

(4,41,0.0),0.0 

S 

0.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

P 

20.00 

0.8294E-04 

0.8294E-04 

0.5604E-05 

0.5604E-05 

S 

20.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

P 

40.00 

0J2282E-02 

0.2282E-02 

0.9586E-05 

0.9586E-05 

S 

40.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

P 

60.00 

0.9161E-04 

0.9161E-04 

0.1403E-04 

0.1403E-04 

S 

60.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

4(  i|(  i|t  41  j|t  s|(  4c  4iifc  %  4t  i|(  %  4*  Xt  %  K(  4c  He  >f(  %  if* 


SIGMA  =  90.00  LAMBDA  =  0.633 

THETA0=  35.00  PfflO  =  0.00 


THETA 

RPl 

RP2 

TPl 

TSl 

THETA 

RSI 

RS2 

TP2 

TS2 

P 

-80.00 

0.1300E-35 

0.1300E-35 

0.4130E-35 

0.4130E-35 

S 

-80.00 

0.1119E-07 

0.1119E-07 

0.1458E-07 

0.1458E-07 

P 

-60.00 

0.1442E-34 

0.1442E-34 

0.2369E-34 

0.2369E-34 

S 

-60.00 

0.1283E-06 

0.1283E-06 

0.1318E-06 

0.1318E-06 

P 

-40.00 

0.8332E-34 

0.8332E-34 

0.844 lE-34 

0.8441E-34 

S 

-40.00 

0.8908E-06 

0.8908E-06 

0.5194E-06 

0.5194E-06 

P 

-20.00 

0.7874E-33 

0.7874E-33 

0.2943E-33 

0.2943E-33 

S 

-20.00 

0.5980E-05 

0.5980E-05 

0.1728E-05 

0.1728E-05 

P 

0.00 

0.7175E-33 

0.7175E-33 

0.7748E-34 

0.7748E-34 

S 

0.00 

0.4272E-04 

0.4272E-04 

0.4612E-05 

0.4612E-05 

P 

20.00 

0.2410E-32 

0.2410E-32 

0.1628E-33 

0.1628E-33 

S 

20.00 

0.2038E-03 

0.2038E-03 

0.8843E-05 

0.8843E-05 

P 

40.00 

0.6630E-31 

0.6630E-31 

0.2785E-33 

0.2785E-33 

S 

40.00 

0.2301E-02 

0.2301E-02 

0.1177E-04 

0.1177E-04 

P 

60.00 

0.2662E-32 

0.2662E-32 

0.4077E-33 

0.4077E-33 

S 

60.00 

0.9495E-04 

0.9495E-04 

0.1164E-04 

0.1164E-04 
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(a)  (b) 

FIGURE  D-1.  Log[ARS]  Versus  Polar  Scattering  Angle  6.  Plotted  is  scattering  in  the  plane  of  incidence 
from  a  bare  surface  (L  =  0)  with  the  incident  beam  of  wavelength  0.633  microns  incident  at  35  degrees.  This 
angle  of  incidence  is  beyond  the  critical  angle  for  specular  transmission.  Positive  and  negative  6  values  are 
for  scattering  azimuth  angles  (|)  of  0  and  180  degrees,  respectively.  The  solid  and  dashed  lines  refer  to  ARS 
when  the  incident  and  scattered  polarization  are  p  -  POL  -*p-  POL  and  s  -  POL  -*s-  POL ,  respectively. 
Reflection  scattering  is  shown  in  (a),  and  transmission  scattering  is  shown  in  (b).  Since  there  is  only  one 
interface,  the  correlated  and  uncorrelated  cases  give  identical  results. 
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Appendix  E 

INPUT  DATA  AND  OUTPUT  RESULTS  FOR  OUT-OF-PLANE  ARS 
FROM  A  FIVE-LAYER  STACK  UPON  A  METALLIC  SUBSTRATE 


OUTPUT  DATA  FILE: 


SIGMA  =  0.00  LAMBDA  =  0.633 

THETA0  =  60.00  PHIO  =  35.00 


P  =»  p-POL  scatt.  S  =>  s-POL  scatt. 

THETA 

RPl 

RP2 

TPl 

TSl 

THETA 

RSI 

RS2 

TP2 

TS2 

P 

-80.00 

0.2290E-06 

0.1626E-05 

O.OOOOE+00 

O.0OOOE+0O 

S 

-80.00 

0.3360E-07 

0.0086E-06 

O.OOOOE+00 

O.0OOOE+0O 

INPUT  DATA  FILE: 

P 

-60.00 

0.6770E-05 

0.1013E-04 

O.OOOOE+00 

O.OOOOE+OO 

S 

-60.00 

0.3337E-06 

0.1760E-05 

O.OOOOE+00 

O.OOOOE+00 

-80.,80.,20. 

P 

-40.00 

0.4994E-06 

0.8653E-05 

O.OOOOE+00 

O.OOOOE+00 

60.,0.633.35.0,5 

S 

-40.00 

0.2639E-08 

0.1664E-05 

O.OOOOE+00 

O.OOOOE+00 

.0005^.0,.00049,0.2 

P 

-20.00 

0.6324E-05 

0.6358E-04 

O.OOOOE+00 

O.OOOOE+00 

(-16.0,0.5).0.0 

S 

-20.00 

0.1478E-05 

0.2828E-04 

O.OOOOE+00 

O.OOOOE+00 

(4.41, 0.0), 0.4 

P 

0.00 

0.3432E-04 

0.5208E-04 

O.OOOOE+00 

O.OOOOE+00 

(2.25,0.0),0.6 

S 

0.00 

0.1683E-04 

0.2553E-04 

O.OOOOE+00 

O.OOOOE+00 

(4.41,0.0),0.4 

P 

20.00 

0.6859E-05 

0.144  lE-03 

O.OOOOE+00 

O.OOOOE+00 

(2.25,0.0),0.6 

S 

20.00 

0.4203E-05 

0.8042E-04 

O.OOOOE+00 

O.OOOOE+00 

(4.41,0.0),0.7 

P 

40.00 

0.1018E-04 

0.4903E-04 

O.OOOOE+00 

O.OOOOE+OO 

(1.0,0.0),0.0 

S 

40.00 

0.1887E-07 

0.1190E-04 

O.OOOOE+00 

O.OOOOE+OO 

P 

60.00 

0.5383E-06 

0.9180E-04 

O.OOOOE+00 

O.OOOOE+00 

S 

60.00 

0.4297E-05 

0.2266E-04 

O.OOOOE+00 

O.OOOOE+OO 

SIGMA  =  90.00  LAMBDA  =  0.633 

THETA0  =  60.00  PHIO  =  35.00 

THETA 

RPl 

RP2 

TPl 

TSl 

THETA 

RSI 

RS2 

TP2 

TS2 

P 

-80.00 

0.4551E-07 

0.2326E-06 

O.OOOOE+00 

O.OOOOE+00 

S 

-80.00 

0.1600E-07 

0.1153E-06 

O.OOOOE+00 

O.OOOOE+00 

P 

-60.00 

0.3337E-06 

0.1760E-05 

O.OOOOE+00 

O.OOOOE+00 

S 

-60.00 

0.1812E-06 

0.1845E-05 

O.OOOOE+00 

O.OOOOE+00 

P 

-40.00 

0.1603E-08 

0.1610E-05 

O.OOOOE+00 

O.OOOOE+00 

S 

-40.00 

0.565  lE-07 

0.1780E-05 

O.OOOOE+00 

O.OOOOE+00 

P 

-20.00 

0.1200E-06 

0.1478E-04 

O.OOOOE+00 

O.OOOOE+00 

S 

-20.00 

0.5996E-06 

0.3106E-04 

O.OOOOE+00 

O.OOOOE+00 

P 

0.00 

0.4013E-05 

0.1217E-04 

O.OOOOE+00 

O.OOOOE+OO 

S 

0.00 

0.8184E-05 

0.2482E-04 

O.OOOOE+00 

O.OOOOE+00 

P 

20.00 

0.3412E-06 

0.4202E-04 

O.OOOOE+00 

O.OOOOE+OO 

s 

20.00 

0.1705E-05 

0.883  lE-04 

O.OOOOE+00 

O.OOOOE+OO 

p 

40.00 

0.1146E-07 

0.1151E-04 

O.OOOOE+OO 

O.OOOOE+OO 

s 

40.00 

0.404  lE-06 

0.1273E-04 

O.OOOOE+OO 

O.OOOOE+OO 

p 

60.00 

0.4297E-05 

0.2266E-04 

O.OOOOE+OO 

O.OOOOE+OO 

s 

60.00 

0.2333E-05 

0.2375E-04 

O.OOOOE+OO 

O.OOOOE+OO 
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e.  DEGREES  B,  DEGREES 


FIGURE  E-l.  Log[ARS]  Versus  Polar  Scattering  Angle  9.  Plotted  is  out-of-plane  scattering  from  a  five-layer 
multilayer  stack  with  the  incident  beam  of  wavelength  0.633  microns  incident  at  60  degrees.  Positive  and 
negative  0  values  are  for  scattering  azimuth  angles  <1)  of  35  and  21 5  degrees,  respectively.  ARS  from  a  correlated 
multilayer  stack  is  shown  in  (a)  where  the  solid  and  dashed  lines  are  for  incident  and  scattered  polarizations 
p -POL  ->  p -POL  and  p -POL  -»  s  -POL,  respectively.  ARS  from  an  uncorrelated  multilayer  stack  is 
shown  in  (b)  where  the  solid  and  dashed  lines  are  for  incident  and  scattered  polarizations  p  -  POL  -^p-  POL 
and  p  -  POL  -^s-  POL ,  respectively.  ARS  from  a  cwrelated  multilayCT  stack  is  shown  in  (c)  where  the  solid 
and  dashed  lines  are  for  incident  and  scattered  polarizations  s  —POL  ->  p -POL  and  s -POL  s -POL, 
respectively.  Finally,  ARS  from  an  uncorrelated  multilayer  stack  is  shown  in  (d)  where  the  solid  and  d^hed 
lines  are  for  incident  and  scattered  polarizations  j  -POL  — >  p  —POL  and  s  —POL  —*  s  —POL,  respectively. 
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Appendix  F 

INPUT  DATA  AND  OUTPUT  RESULTS  FOR  ARS  FROM  A  SEMI-OPAQUE 
THIN  METALLIC  FILM  INCLUDING  SURFACE  PLASMA  RESONANCE 


OUTPUT  DATA  FILE: 


SIGMA  = 

0.00  LAMBDA  = 

0.633 

THETA0  =  29.60  PfflO  = 

0.00 

P  =»  p-POL  scatt.  S  =>  s-POL  scatt. 

THETA 

RPl 

RP2 

TPl 

TSl 

THETA 

RSI 

RS2 

TP2 

TS2 

P 

-80.00 

0.9282E-07 

0.4515E-06 

0.7616E-05 

0.1086E-04 

S 

-80.00 

0.9677E-37 

0.7410E-35 

0.2524E-34 

0.5526E-34 

INPUT  DATA  FILE: 

P 

-60.00 

0.1013E-05 

0.4664E-05 

0.2858E-04 

0.4199E-04 

S 

-60.00 

0.5628E-36 

0.1005E-33 

0.3041E-33 

0.6638E-33 

-80.,80..20. 

P 

-40.00 

0.1272E-04 

0.4540E-04 

0.5558E-04 

0.8714E-04 

29.6,0.633,0.0,1 

S 

-40.00 

0.1164E-35 

0.1076E-32 

0.1438E-32 

0.3125E-32 

.0005,2.0,.00049,0.2 

P 

-20.00 

0.2593E-04 

0.1451E-03 

0.9204E-04 

0.1620E-03 

(1.0,0.0),0.0 

S 

-20.00 

0.6472E-33 

0.1907E-31 

0.5064E-32 

0.1096E-31 

(-16.0,0.5),0.025 

P 

0.00 

0.5561E-04 

0.1221E-02 

0.1103E-03 

0.2383E-03 

(4.41,0.0),0.0 

S 

0.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

P 

20.00 

0.6053E-04 

0.5371E-02 

0.7901E-04 

0.2540E-03 

S 

20.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

P 

40.00 

0.3037E-04 

0.6166E-02 

0.3026E-04 

0.2375E-03 

S 

40.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

P 

60.00 

0.1903E-04 

0.1190E-02 

0.4425E-05 

0.5173E-03 

S 

60.00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

O.OOOOE+00 

SIGMA  =  90.00  LAMBDA  = 

0.633 

THETA0  =  29.60  PHIO  = 

0.00 

THETA 

RPl 

RP2 

TPl 

TSl 

THETA 

RSI 

RS2 

TP2 

TS2 

P 

-80.00 

0.9427E-35 

0.2177E-34 

0.1645E-33 

0.3375E-33 

S 

-80.00 

0.1480E-07 

02585E-07 

0.3452E-08 

0.4845E-07 

P 

-60.00 

0.9127E-34 

0.2153E-33 

0.6079E-33 

0.1315E-32 

S 

-60.00 

0.1838E-06 

0.3349E-06 

0.4370E-07 

0.5827E-06 

P 

^0.00 

0.7838E-33 

0.1847E-32 

0.1138E-32 

0.2777E-32 

S 

-40.00 

0.161  lE-05 

0.3236E-05 

0.2219E-06 

0.2748E-05 

P 

-20.00 

0.5873E-32 

0.7613E-32 

0.1738E-32 

0.533  lE-32 

S 

-20.00 

0.2058E-04 

0.4852E-04 

0.8282E-06 

0.9656E-05 

P 

0.00 

0.1616E-32 

0.3548E-31 

0.3205E-32 

0.6925E-32 

s 

0.00 

0.1545E-03 

0.3602E-03 

0.2143E-05 

0.2442E-04 

p 

20.00 

0.1759E-32 

0.1561E-30 

0.2296E-32 

0.738  lE-32 

s 

20.00 

0.6590E-03 

0.1554E-02 

0.3543E-05 

0.4131E-04 

p 

40.00 

0.8825E-33 

0.179  lE-30 

0.879  lE-33 

0.6900E-32 

s 

40.00 

0.5025E-03 

0.1009E-02 

0.4206E-05 

0.5210E-04 

p 

60.00 

0.5530E-33 

0.3456E-31 

0.1286E-33 

0.1503E-31 

s 

60.00 

0.7567E-04 

0.1379E-03 

0.7709E-05 

0.1028E-03 
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FIGURE  F-1.  Log[ARS]  Versus  Polar  Scattering  Angle  0.  Plotted  is  scattering  in  the  plane  of  incidence  for 
a  single-layer  case  (L  =  1)  with  the  incident  beam  of  wavelength  0.633  microns  incident  at  29.6  degrees. 
Positive  and  negative  0  vadues  are  for  scattering  azimuth  angles  0  of  0  and  180  degrees,  respectively.  The 
solid  and  dashed  lines  refer  to  ARS  when  the  multilayer  stack  interface  roughness  is  correlated  and  uncorrelated, 
respectively.  Reflection  scattering  is  shown  in  (a)  and  (b),  and  transmission  scattering  is  shown  in  (c)  and 
(d).  The  incident  and  scattered  polarization  for  these  four  plots  are  (a)  p-POL p-POL,  (b) 
j -POL  -4  5 -POL,  (c)  p -POL  p -POL,  and  (d)  s -POL  s -POL. 


4. 
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